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Absolute parameters of AE For — a highly active detached 
binary of late K type* 
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ABSTRACT 

We present photometric and spectroscopic analysis of AE For - a detached eclipsing 
binary composed of two late K dwarfs. The masses of the components are found 
to be 0.6314±0.0035 and 0.6197±0.0034 M Q and the radii to be 0.67±0.03 and 
0.63±0.03 Rq for primary and secondary component, respectively. Both components 
are significantly oversized compared to theoretical models, which we attribute to their 
high activity. They show H a , H^, H 7 , Ca H and Ca K in emission, and are heavily 
spotted, causing large variations of the light curve. 
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1 INTRODUCTION 

The 10th magnitude st ar AE For (CD -25 1273; HIP 14568) 
was first classified by IStephensonl ll 19861) as a K4 dwarf. 
In the Hipparcos catalogue ( ESA] 1 19971 ) it was listed as 
a new eclipsing variable of Algol type with a period of 
0.918235(8) d, and a parallax of 32.10±1.78 mas. The latter 
was revised in 2007 to the presently adopted 31.8±1.96 mas 
Jvan Leeuw cn 2007) , corre s pondi ng to 31.5±1.9 pc. 

(jlZIS, Reid fc HawlevI (|2002T l found the star to be a 
double-lined binary with Balmer lines in emission. Emis- 
sion in Balmer an d Ca II H & K lines was also observed by 
iGrav et al.l l|2006h . who estimated the spectral type of AE 
For at K9 Ve, and classified it as a very active system. A 
slightl y earlier spectra l type (K7 Ve) was assigned to AE For 
by I Torres et al.l (|2006T l. who also noted its high activity. Ac- 
cording to Zasche, Svoboda & Uhlaf (2012; hereafter ZSU), 
the system contains a brown dwarf with a minimal mass of 
~47 Mj up on an eccentric orbit of ~7 years. 

In the solar neighborhood about 7 per cent of de- 
tached eclipsing binaries are X-ray emitters 
2008]! . With its L x = 9.4 x 10 29 erg 
199a) and F x /F opt 
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* Based in part on data obtained at the Las Campanas Observa- 
tory. 
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is a prominent memb er of this group. In the catalogue of 
ISzczygiel etHI (|2008l l it is one of the closest X-ray active 
binar ies, whose Lx/L opt ratio reaches 0.01 (|Fischer et al.l 
ll999|V An X-ray flaring act i vity o f the system was reported 
by iFuhrmeister fc Schmittl (|2003T l. It is also a bright IR 
source, with 2MASS magnitudes J = 7.52, H = 6.85 and 
K = 6.66. 

For about a decade it has been known that active dwarfs 
of K and M type tend to be larger and cooler than the the- 
ory predicts. While possible solutions of this problem have 
been proposed (see e.g. Morales, Ribas & Jordi 2008 and 
references therein), it is certainly worthwhile to enlarge the 
relevant observational databasis. The best opportunity for 
this is offered by detached binaries on the lower main se- 
quence, which allow to determine masses and radii of their 
components with an accuracy better than 1 per cent. AE 
For is clearly one of such systems, however until very re- 
cently neither the light curve nor the velocity curve of this 
potentially interesting binary has been studied. A prelim- 
inary light-curve solution has been derived by ZSU, who 
concluded that a spectroscopic analysis was needed to con- 
firm the physical parameters of the components to a higher 
accuracy. 

In the present paper we obtain and analyze the veloc- 
ity curve of AE For, and refine the photometric solution 
of ZSU based on additional observations. The photometric 
and spectroscopic data are described in Sects. 2 and 3. The 
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Table 1. List of light curves used in this paper for photometric 
solutions of AE For. 



I s * day 


last day 


Filter 


Number of 


Ref. 


[HJD-2400000] 




data points 




55174 


55184 


Johnson B 


160 


1 


55174 


55184 


Johnson V 


482 


1 


55568 


55578 


Johnson B 


173 


2 


55568 


55578 


Johnson V 


175 


2 


55568 


55578 


Cousins R 


173 


2 



1: our own data; 2: IZasche et al. 

analysis of the data is detailed in Sect. 4, and its results are 
discussed in Sect. 5. 




2 PHOTOMETRIC OBSERVATIONS 

The observational material consists of our own BV data 
and BVR data of ZSU (see Table [T] for details). Our data 
were collected during 11 nights between 2009 Dec 8/9 and 
Dec 18/19 with the 1-m Swope telescope at Las Campanas 
Observatory (LCO), Chile. All the nights were clear, with 
the seeing between 1.1 and 2.5 arcsec. AE For was monitored 
with the 2048 x 3150 pixel SITe3 CCD camera at a scale of 
0.435 arcsec/pixel. We collected 482 frames through the V 
filter with exposure times 5-10 s (depending on the seeing), 
and 160 frames in the B filter with exposure times 8-25 s. 
All images were de-biased and flat-fielded within the 

For photo metric measurements the Daophot package 
(Stetson 1987) was employed, and differential aperture 
photometry was extracted using two stars in the vicin- 
ity of AE For. On every night AE For was observed we 
also observed NGC 2204. For six red clump stars of the 
latter we found the offset between instrumental magni- 
tudes and the stand ard BV magnitudes determined by 
I Rozyczka et al.l (|2007T ) . Using that offset together with the 
average extinction for LCO (ky = 0.14 and ks = 0.24; see 
iMininti, Claria fc Gomez! 1 19891 ). and accounting for a 0.04 
difference in air-mass between AE For and NGC 2204, we 
transformed the instrumental light curves of the variable to 
the standard BV system. Color terms of the transformation 
were accounted for. We estimate the error of this transfor- 
mation at ±0.02 mag, its main source being the uncertainty 
of the zero-point for NGC 2204. The formal error of the 
differential photometry is about ten times smaller. 

The last three rows in Table [1] refer to 
BVR measurements performed by ZSU in Jan- 
uary 2011 at the South African Astronomical 
Observatory, (their data are accessible online at 
|http : / / viz ier . cf a .harvard . edu / viz- bin / VizieR?-sourcc= J/A+A/53 
Accessible online were also Hipparcos (ESA 1997), AS AS 
l|Poimanskil 120021 ) . and Pi of the Sky (|Burd et all 120051 ) 
light curves. However, their quality was too poor to use 
them for photometric solutions. 

Based on almost 70 primary and secondary minima 



Figure 1. Light curves of AE For obtained in December 2009 
at Las Campanas Observatory. To save the plotting space the 
B-curve is shifted upward by one magnitude. 



observed by various authors, ZSU derived the following 
ephemeris of AE For: 

tfJAnin = 2452605.97070(35) + 0.91820943(12) x E (1) 

The final light curves obtained from the data collected at 
LCO and phased with this ephemeris are shown in Fig. [T] 
Evident is a strong asymmetry, which we attribute to the 
activity of the system, manifesting itself by large spot(s) on 
at least one of the components (see also Sects. [3] and 3}. 
For the short flat part of the U-curve which begins right 
after the secondary eclipse (see Fig. [TJ we obtained V = 
10.27±0.02 mag. To within the errors, our (B — V) index 
was constant throughout the orbital period, and equal to 
1.35±0.03 mag. For the same part of the light curve the 
online data of ZSU yield V = 10.23±0.01 mag and (B-V) = 
1.34±0.02 mag. Marginal differences between our results and 
theirs may be due to the spot-related variability of at least 
one component of AE For. 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the AURA, Inc., under coopera- 
tive agreement with the NSF. 



3 SPECTROSCOPIC OBSERVATIONS 

The spectroscopic data were collected with the Echelle Spec- 
trograph on the 2.5-m Irenee du Pont telescope at LCO. The 
spectrograph, equipped with the SITe2K CCD camera, was 
working at a slit width of 1.0 arcsec, providing a resolution of 
^45,000. The observations were performed in two runs: four 
nights from 2009 Nov 23/24 to Nov 26/27 and five nights 
from 2010 Dec 19/20 to Dec 23/24. All the nights were clear 
but a half of one night with thin cirrus clouds. For most 
fj^WoW) 1 cent) exposures the airmass was smaller than 1.4. 
uuring the observations pairs of 480-s exposures of the sci- 
entific target were made, interlaced with a 90-s exposure of 
a ThAr lamp spectrum. 

The observations were reduced within the IRAF 
ECHELLE package. After bias and flat-field correction, each 
pair of the frames was combined into a single frame, allowing 
for the rejection of cosmic ray hits. Altogether, 31 reduced 
spectra were obtained, extending from ~4000 to ~7000 A. 
All of them showed strong Balmer emission lines, and in 
nearly all of them the Nal doublet at 5890 and 5896 A was 
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Table 2. Barycentric radial velocities of AE For phased according 
to the ephemeris given by equation pj. 



HJD-2455000 


Vp [km s 1 ] 


v s [km s x ] 


phase 


1 ^791 3 
ioy .O / ZlO 


y.oo 


110 40 


n OfifiAA 


ioy .ooooz 


98 Q1 

-zo. yi 


1 A8 Q^i 

i^to.yo 


U. lOooO 


1 7077^ 

ioy . iU( iu 


^9 
-oz.yo 


1 77 QA 
1 I t .y4: 


21414 


ioy . i / ooy 


-oo.oo 


1 7^ (17 


n 98fi1 7 
U.ZOOl / 


1 7RA1 

ioy . i o^txo 


^1 ^8 

-Jl. OO 


1 7*3 1 Q 

i i o.iy 


n 9Q7^fi 
U.zy / ou 


lOU.OZOOU 


1 A A7 


1 LA R1 


n i n^fin 

U. 1UODU 


luu.ozyoo 


^3 Q8 

-oo.yo 


1 7^ 8*3 
i / o.oo 


21841 


1 fin flAfiQQ 


^fi 77 
-OD. / / 


1 7Q 90 

± / y.zu 


n 9^7nA 
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-O 1 .10 


1 7Q ^8 
1 i y.oo 


n 9fi ,: t89 


1 fiO 7^^QQ 

iou. i ooyy 


'in ^8 


1 1^1 Q7 

ioi. y i 


U.OOO i 


1U1 .Ul I DO 


^9 8n 
-oz.ou 


174 11 


U.zyoo / 




1 17 
o. 1 / 


"1 "1 7 71 
1 1 / . ( 1 


n A9i ^ 

U.'tZloO 


1 fi9 ^'39 7R 

lOZ.OOZ / O 


^3 ^7 


1 7 A 3 A 


n 9Q081 
u.zyuoi 


1 R9 ^AQ&fi 


aq nn 
-4ty .uu 


1 I u.uu 


n ^nQ9n 
u.ouyzu 


1 (\0 (^07^7 
1UZ.UU I O I 


-24 52 


1 A^ fi^ 


n ^7997 
U.O 1 ZZ / 


1 fi9 fiA1 9Q 

loz.o^izy 


-o.oy 


1 9fi fi'? 


n An&QQ 
u.^tuoyy 


1 fi9 R1 9A1 
1UZ . OlZ^il 


126 34 


7 AQ 


n AnAfi^ 


162.82922 


138.43 


-20.31 


-0.38634 


550.55538 


143.22 


-24.08 


-0.12303 


550.57283 


129.29 


-10.95 


-0.10402 


550.58045 


124.95 


-6.27 


-0.09573 


551.73601 


-38.74 


161.16 


0.16277 


551.75427 


-47.05 


168.90 


0.18265 


552.68821 


-50.75 


173.13 


0.19979 


552.70336 


-54.04 


175.98 


0.21629 


552.71858 


-56.71 


178.77 


0.23286 


553.59737 


-49.76 


170.50 


0.18993 


553.67500 


-54.54 


177.94 


0.27448 


553.73276 


-39.80 


162.90 


0.33738 


554.61505 


-50.65 


172.56 


0.29827 


554.67243 


-27.41 


150.57 


0.36076 



blended into a very broad (~15 A, or ~760 km s _1 ) absorp- 
tion feature. 

Radial velocities were measured using an implementa - 
tion of the b roadening func ti on fo rmalism l|Rucinskil 12003 ) 
described by iKaluznv et all (|2006h FI Since including emis- 
sion lines or broad absorption features would decrease the 
accuracy of velocity measurements, a wavelength range de- 
void of such features had to be selected. Shortward of H 7 
the spectra were too noisy, so that our choice was reduced 
to four segments limited by H 7 H^, Nal doublet, H a and O2 
b-band. We chose the segment between H/3 and Nal dou- 
blet, extending from 4870 to 5845 A. It was the longest one 
of the four, and the only one in which the mean S/N ratio 
was larger than 20 for all the spectra. The solar-scaled syn- 
thetic sp ectrum for T = 400 K and logg = 5.0 from the 
library of lCoelho et~all |2005l ) served as a template. 

The observed velocity curve was fitted with a nonlin- 
ear least-squares solution, using a spectroscopic data solver 
written and kindly provided by Guillermo Torres. Because to 
within observational errors secondary minima are separated 
from the primary ones by half the period, a circular orbit 
was assumed, and the eccentricity e was fixed at (the same 
assumption was adopted by ZSU). The resulting barycentric 

2 The software package used in this paper is freely accessible at 
http://uscrs.camk.cdu.pl/pych/BF/ 
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Figure 2. Velocity curve of AE For based on the data contained 
in Table [2] The rms residual velocity calculated from all (i.e. 
primary's and secondary's) residua shown in the bottom panel is 
equal to 1.13 km s -1 . Phase corresponds to the center of the 
primary photometric minimum. 

radial velocities are listed in Table [2] and plotted in the up- 
per panel of Fig. [2] together with the fitted velocity curve. 
The mean error of velocity measurement, estimated from 
the residual velocities shown in the bottom panel of Fig. [2] 
is ±1.13 km s _1 . Its rather large value may result from line 
blending and/or asymmetries (the latter being caused by 
large spots; see also Sects. [2] and UJ. We repeated the mea- 
surements using a synthetic spectrum for T = 4250 K to 
find that the changes in the observed velocities were much 
smaller than the scatter of points around the fit. We also 
tried out the remaining three wavelength ranges, however in 
all cases the fit had a lower quality (i.e. the residual veloci- 
ties were larger). 

The derived orbital parameters are listed in Table [3] 
together with formal 1-cr errors returned by the fitting rou- 
tine. Note that, despite the rather large residual velocities, 
the solution is fairly accurate - relative errors in M sin 3 i and 
Asini amount to 0.5 and 0.15 per cent, respectively. 



4 LIGHT CURVE ANALYSIS AND SYSTEM 
PARAMETERS 

The analysis of the light c urves was performed with the 
PHOEBE implementation (|Prsa fc Zwitterl l2005h of the 
Wilson-Devin ney (WD) model l|Wilson fc Devinnevl Il97ll ; 
IWilsonl Il979l ). The PHOEBE/ WD package utilizes the 
Roche geometry to approximate the shapes of the stars, 
uses Kurucz model atmospheres, treats reflection effects in 
detail, and, most importantly, allows for the simultaneous 
analysis of B and V data. Appropriately for stars with con- 
vective envelopes, we adopted gravity darkening coefficients 
9p = 9a = 0.32 and bolometric albedos A p = A s = 0.5. 
The effects of reflection were incl uded, and a l ogarit hmic 
limb-darkening based on tables by Ivan Hammd (Il993l ) was 
used as implemented in PHOEBE 031a. Following ZSU, full 
synchronization of both components was assumed. 

ZSU fitted the lig ht curves of A E For with the help of 
the program ROCHE jPribuUdfedolh . Keeping in mind that 
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Table 3. Orbital parameters 



Parameter 


Unit 


Value 


7 


km s _1 


59.95±0.23 


K p 


km s~ 1 


116.92±0.33 


K s 


km s _1 


119.12±0.33 


e 




0.0 a 


(Tp 


km s _1 


1.08 




km s _1 


1.18 


Derived quantities: 






A sin i 


Re 


4.2820±0.0069 


M p sin 3 i 


M 


0.6314±0.0034 


M s sin 3 i 


M© 


0.6197±0.0033 



a Assumcd in fit 



their ROCHE and PHOEBE fits of GK Boo were markedly 
different (although not at a statistically significant level), we 
decided to check if our photometric solution of AE For will 
be compatible with theirs. We fed their photometric data 
into PHOEBE, and iterated upon inclination i, secondary's 
temperature T s , primary's radius R p and secondary's radius 
R s , using values obtained by ZSU as initial values. The mass 
ratio q was fixed at a value of 0.98522 calculated from Ta- 
ble[3] and the orbital separation A - at 4.28996 Rq obtained 
from Asini in Table for the inclination i — 85.51° found 
by ZSU. Following ZSU, we kept the temperature of the pri- 
mary T p fixed at 4100 K. To within the errors, this value 
agrees with 4 055±71 K resulting fr o m the T e // — (B — V) 
calibration of I Ramirez fc Melende j l|2005l ). and marginally 
agrees wit h 4005±92 K resulting from the analogous cal- 
ibration of IWorthev &: Lee] (|20 1 lT ) - The latter two values, 
however, should be treated with some caution, since AE For 
emits most energy in the near-IR range beyond B and V 
passbands, at 2MASS magnitudes J = 7.52, H = 6.85 and 
K — 6.66. An infrared calibration would be more appro- 
priate; unfortunately we have not found any such relation 
extending up to (J — K) = 0.86. The problem of the tem- 
perature is further discussed in Sect. [5] 

The results of fitting are listed in columns 3 and 4 of 
Table [4] Column 3 contains the original values obtained by 
ZSU together with their errors. The corresponding synthetic 
light curves produced by PHOEBE without any iterations 
turned out to be fairly well fitting, with a standard devia- 
tion of V-residuals equal to 17 mmag (similar values were 
obtained for B and R bands). The iterations improved the 
quality of the fit (ov was reduced to 13 mmag), but the 
parameters did not change in a statistically significant way 
(column 4). We performed an additional iteration with the 
following spots introduced on both components at a latitude 
of 90° (i.e. at the equator in PHOEBE's convention): one on 
the primary at a longitude of 180° with a radius of 20° and 
temperature factor of 0.98; two on the secondary at longi- 
tudes = 125° and 330° , both of them with a radius of 30° 
and temperature factor of 1.02. The residuals became more 
symmetric (see Fig. [3| and av was reduced to 12 mmag, but 
no statistically significant corrections to system parameters 
were obtained (column 5 of Table[4]). We note that by "spot" 
we mean a fairly large region on the surface of a star, whose 
mean temperature is elevated or reduced due to an excess 
of genuine stellar spots with much smaller sizes. 

For the assumed M s = 0.5A/©, ZSU derived R s = 



Table 4. Photometric parameters derived from ZSU data. The 
lower bounds for component masses derived from the orbital so- 
lutions are 0.6314±0.0034 and 0.6197±0.0033 Mq, respectively, 
for the primary and the secondary. 



Parameter 


Unit 


Original 


Our fit 


Our fit 


Dark 






ZSU 




with spots 


pole 


1 


2 


3 


1 


5 


6 


i 


deg 


86.51(31) 


86.84 


86.71 


85.47 


T s 


K 


4065(48) 


4083 


4052 


4041 


Rp 


Rq 


0.66(10) 


0.71 


0.71 


0.661 


Rs 


R 3 


0.52(8) 


0.53 


0.52 


0.605 


(L p )v 


% 


63.1(1.2) 


65.4 


66.8 


57.2 


(L p )b 


% 


63.2(1.3) 


65.6 


67.3 


57.9 


a v 


mmag 


17 


13 


12 


13 
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Figure 3. Residuals of fits listed in columns 3—5 of Table [4] 
Top: original ZSU: middle: PHOEBE-iterated starting from ZSU 
parameters, bottom: same as middle, but with spots added on 
both components (see text for details). 

0.52i?Q, and our PHOEBE iterations confirmed this result. 
However, while a radius of 0.52_R© may be compatible with 
the mass of M s = 0.5M© assumed by ZSU, it is excluded 
by our spectroscopic solution, in which M s > 0.6M©. This 
is because within a broad range of metallicity (—0.5 ^ 
[Fe/H] 0.5) and a-enhancement (-0.2 < [a/Fe] sC 0.4) 
a main-se quence star w ith M=O.6M0 never becomes that 
small (see Po tter et al.ll2008l . and the DSED evolutionary 
tracks at http://stellar.dartmouth.edu/models/index.htm). 
Another argument in favor of a larger secondary is based on 
rotational velocity fits which are performed automatically 
within the BF formalism during radial velocity measure- 
ments. Upon averaging velocities fitted in 2009 and 2010 we 
got Urot.s = 37.46±2.40 km s -1 for the primary and v Iot ,s = 
37.07±2.20 km s _1 for the secondary. For the assumed syn- 
chronous rotation these values imply R p = 0.680 ± O.O44i?0 
and R s = 0.673 ± 0.0407?©. 

The simplest way to remove this contradiction is to 
force a smaller inclination i: one may expect that Rs will 
then increase, so that depths and widths of the minima are 
preserved. This may be achieved, for example, by placing 
a dark spot around the uneclipsed pole of the primary. We 
performed an experiment, in which the circumpolar region of 
the primary with a radius of 50° had a temperature factor of 
0.95. As shown in column 6 of TablelH the iterations yielded 
a solution with i smaller by more than one degree and R s 
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Table 5. Photometric parameters derived from our data 



Parameter 


Unit 


PO 


P0 


Plh 


Pic 


Mean 






SO 


Slh2c 


S2c 


Slhlc 


col. 4-6 


1 


2 


3 


4 


5 


6 


7 


i 


deg 


85.72 


85.75 


85.51 


85.43 


85.6±0.2 


T a 


K 


4016 


4055 


4061 


4050 


4055±6 


Rp 


Hq 


0.705 


0.694 


0.665 


0.650 


0.67±0.03 


Rs 


R : 


0.603 


0.607 


0.640 


0.654 


0.63±0.03 


(L P )v 


% 


61.6 


58.8 


53.8 


52.1 


54.9±4.4 


(L p )b 


% 


62.5 


59.4 


54.2 


52.7 


55.4±4.4 


(Tv- 


mmag 


20.1 


7.0 


6.5 


6.5 
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Figure 4. Residuals of fits listed in columns 3 and 4 of Table [5] 
Top: model P0S0 (without spots); bottom: model P0Slh2c (three 
spots on the secondary; see Table [5] for details). 



Table 6. Parameters of spots for models from Table [5] 



Model 


component 


latitude 


radius 


temp, factor 


P0Slh2c 


S 


60 


40 


0.99 




s 


100 


30 


1.02 




s 


230 


30 


0.95 


PlhS2c 


P 


280 


30 


1.02 




s 


60 


10 


0.99 




s 


230 


30 


0.95 


PlcSlhlc 


P 


240 


40 


0.99 




s 


100 


30 


1.02 




s 


230 


30 


0.95 



larger by almost 0.1 Rq compared to the model without a 
darkened pole. 

At this point we faced the problem whether to include 
clearly unphysical solutions in final estimates of system pa- 
rameters, and we decided to reject them. Apparently, ZSU 
collected their data when the system was in a particular 
state making the analysis very uncertain, and to remain on 
the safe side we calculated the parameters of AE For based 
on our photometry only. 

Photometric solutions based on our own data are listed 
in Table [5] Four models are shown, identified by the num- 
bers of hot (h) and cold (c) spots on the primary (P) and 
the secondary (S). In each spotted model there are three 
spots centered at a latitude of 90°. The remaining parame- 
ters of the spots are listed in Table [6] The fit without spots 
(P0S0) is poor, as indicated by a large ov ■ Its low quality 
is clearly seen in Fig. 2] in which the residuals show sys- 
tematic deviations with an amplitude exceeding 40 mmag. 
Introducing spots makes the fit nearly ideal. Not surpris- 
ingly, however, Table [S] shows that different spot arrange- 
ments result in different system parameters. The encourag- 
ing finding is that for all models (even for the unspotted 
one) the inclination angle is by about one degree smaller 
than that obtained from ZSU data, and the radius of the 
secondary is consistently larger than 0.6Rq. The last col- 
umn of Table [5] contains averaged values of spotted model 
parameters together with formally calculated errors. Since 
the distribution of spot-dependent model parameters is most 
probably non-Gaussian, these errors should be treated with 
caution, and we give them only to approximately illustrate 
the uncertainties resulting from the nonexistence of a unique 
configuration of spots. 

Independently of their spot-related uncertainty, the 
models have Gaussian errors related to nonzero residuals, 



which we estimated using a Monte Carlo procedure written 
in PHOEBE-scripter. Briefly, the procedure replaces the ob- 
served light curves B a i, s and V b 3 with the fitted ones Bf 
and Vf, generates 20000 Gaussian perturbations SBf and 
SVf such that the standard deviation of each of them is 
equal to the standard deviation of the residuals Bf — B b s 
or Vf — Vots, and for each perturbation performs PHOEBE 
iterations on Vf + SVf and Bf + SBf. In all cases the Gaus- 
sian errors turned out to be much smaller than the formal 
errors of averaged parameters given in column 7 of Table [5] 
The final absolute parameters of AE For are given in 
Table [7] (bolometric magnitudes in the last two rows are av- 
erages of the values taken directly from PHOEBE output). 
The errors of A, M p and M s include the inclination uncer- 
tainty. As we explained above, the errors of the remaining 
parameters are but an approximate illustration of the un- 
certainties related to the presence of spots. 



5 DISCUSSION AND CONCLUSIONS 

The fits described in Sect. [3] indicate that an unavoidable 
consequence of spot activity is the degradation of the photo- 
metric solution: based on the available data, the radii of the 
components cannot be calculated with an accuracy better 
than ~5 per cent. As we have shown in Sect. 2] the asym- 
metric and highly variable light curve of AE For can even 
generate nonphysical solutions in which one of the stars is 
unrealistically small. The only way to improve the accuracy 
of system parameters is via collecting more photometric data 
from various seasons and careful fitting of each light curve 
with various arrangements of spots (see e.g. iRozvczka et alj 
2009). Alternatively, many light curves from various seasons 
might be averaged in hope that the mean curve would be 
free from spot-caused irregularities. 

Another potential source of uncertainties in our final pa- 
rameters is the effective temperature of the primary, which 
we fixed at 4100 K (i.e. a value adopted by ZSU). Luck- 
ily, the temperature only weakly affects i, R a and R p - we 
checked that they all change by less than 0.2 per cent for a 
±150 K change in T p . It is thus possible to compare the de- 
rived radii with stellar models, and such a comparison with 
the synthetic main sequence represented by the 1 Gyr Dart- 
mouth isochrone for [Fe/H] = 0.0 and [a/Fe] = 0.0 is shown 
in Fig. [5] Each component of AE For is clearly oversized 
vs. a main-sequence star of the same mass. This conclusion 
holds for a broad range of chemical composition parameters 
(-0.5 < [Fe/H] < 0.5; -0.2 < [a/Fe] < 0.4). 
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Table 7. Absolute parameters of AE For 



Parameter 


Unit 


Value 


A 


R, 3 


4.2820±0.0069 


i 


deg 


85.6 ±0.2 a 


e 




b 


M p 


Mq 


0.6314±0.0035 


M s 


Mq 


0.6197±0.0034 


Rp 


Rq 


0.67±0.03 a 


R s 


Rq 


0.63±0.03 a 


T P 


K 


4100 b 


T s 


K 


4055±6 a 


M^° l 


mag 


7.15±0.09 a 


M lol 


mag 


7.32±0.11 a 



a Approximate errors due to spots (see text) 
b Assumed in fit 



VI [VI. j 

Figure 6. Location of AE For components in the (M — My) 
plane. The magnitudes are derived from the observational data. 
Solid line: the same isochrone as in Fig. [5] 




M [MJ 

Figure 5. Location of our final solution in the [M — K) plane. 
Solid line: main sequence represented by the 1 Gyr Dartmouth 
isochrone for solar abundances. 



While active components of eclipsing binaries appear to 
be larger and cooler than inactive sin gle stars of the sam e 



mass, they have a similar luminosity (|Morales et ai1 [2008L 
To check if this holds for our binary, we calculated absolute 
magnitudes of the components in V-band from the observed 
magnitude and parallax of the system, using the contribu- 
tion of the primary to the total light from Table [7] The 
location of the components in the [M — Mv) plane is shown 
in Fig. [5] together with the main sequence represented as 
before by the 1 Gyr Dartmouth isochrone for solar abun- 
dances. As one can see, to within the errors both stars do 
indeed belong to the main sequence. In the following we as- 
sume that they have main-sequence luminosities, and based 
on that assumption we estimate their temperatures. 

Fig. indicates that the solution with T p = 4100 K 
is significantly too luminous. Only at T p = 3900 K (for 
which T a — 3860 K) do both stars align with the main 
sequence, suggesting that the temperature assumed by 
Z SU and the temperat u res d erive d from the calib r ations 
of iRamirez fc Melendej (|2005h and IWorthev fc Led (|201ll ) 
are overestimated. Of course, this finding should b e verified, 
prefe rably by using disentangling software (see e.g. lHadraval 
l2009h . Further spectroscopic observations are needed for 
that purpose, collected during low-activity periods. When 
the system is active, all temperature estimates, whether 
based on spectral fitting or photometric calibrations, are 




VI [Mj 

Figure 7. Location of our solutions in the (M — L) plane. Models 
with T p = 4100 and T P = 3900 K are identified by filled and open 
circles, respectively. Solid line: the same isochrone as in Fig. [5] 



likely to be flawed because of large spotted areas present 
on both components. 

We have thus proved the consistency of the assumption 
that the components of AE For are main-sequence stars. 
However, given the large errors of My, a possibility that 
they have not yet settled on the main sequence should also 
be explored. If they were indeed in the contraction phase, 
their large sizes could be at least partly caused by evolution- 
ary effects. Such a possibility is indicated by t he Li 6708 A 
line id entified in the spectrum of AE For by iTorres et alj 
(2006). With an equivalent width of 80 mA it falls right 
in the middle of the range observed for Pleiades members 
with the same (V — I), suggesting that the system is rather 
yo ung. Unfortu n ately, our efforts to repeat the measurement 
of lTorres et all ( 2006) were unsuccessful. We do not claim 
that the line is not there: our spectra may be too noisy, or 
the broadening is too strong, or both. It is w orth mention - 
ing, however, that based on kinematic criteria, lE gg enl dl99fi 
assigned AE For to the Hyades superclust er whose age is 
estimated at ~0.6 Gyr (|Montes et al.l l200lL If this assign- 
ment is correct, then both components of the binary must 
have ended their pre- main sequence evolutionary phase long 
ago. Again, high-quality spectra, collected preferably dur- 
ing eclipses, are needed to resolve the Lithium (and age) 
problem. 
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Yet another explanation for the oversized components 
involves the third body discovered by ZSU. In a triple sys- 
tem consisting of a binary orbited by a distant companion 
both the eccentricity of the binary and the inclination of 
the third body's orbit execute periodic Kozai oscillations. 
If at some phase of the oscillation cycle eb becomes suffi- 
ciently large for tidal friction to dissipate the orbital en- 
ergy of the bi nary, then the binary gradu ally tightens its 
orbit (see e.g. iFabrvckv fe Tremaind I2007I . and references 
therein). The associated tidal heating could temporarily in- 
flate the components; however, since nearly all close bi naries 
are members of triple systems (jTokovinin et al.l [2006). and 
the infla tion effect is limite d to low-mass stars of K and 
M type l|Morales et alj|20 08l, we consider this explanation 
much less likely than the remaining two (especially than that 
related to stellar activity). 

The above discussion indicates that AE For is a truly 
interesting object which clearly deserves closer attention. A 
systematic study of this system would certainly bring valu- 
able information concerning the activity of late-type stars. 
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Table 8. A sample of the B-lightcurve of AE For. 



HJD-240000 B [mag] 



55174.610258 11.629 

55174.629343 11.680 

55174.650730 11.994 

55174.673195 11.956 

55174.695058 11.646 

55174.715983 11.619 



Table 9. A sample of the V-lightcurve of AE For. 



HJD-240000 


V [mag] 


55174.662443 


10.745 


55175.581451 


10.744 


55174.664191 


10.731 


55175.583256 


10.733 


55174.665765 


10.725 


55174.667675 


10.697 



ONLINE MATERIAL 



